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1 Overview 



The planar dimer model is, from one point of view, a statistical mechan- 
ical model of random 2-dimensional interfaces in M 3 . In a concrete sense 
it is a natural generalization of the simple random walk on Z . While the 
simple random walk and its scaling limit, Brownian motion, permeate all of 
probability theory and many other parts of mathematics, higher dimensional 
models like the dimer model are much less used or understood. Only re- 
cently have tools been developed for gaining a mathematical understanding 
of two-dimensional random fields. The dimer model is at the moment the 
most successful of these two-dimensional theories. What is remarkable is 
that the objects underlying the simple random walk: the Laplacian, Green's 
function, and Gaussian measure, are also the fundamental tools used in the 
dimer model. On the other hand the study of the dimer model actually uses 
tools from many other areas of mathematics: we'll see a little bit of algebraic 
geometry, PDEs, analysis, and ergodic theory, at least. 

Our goal in these notes is to study the planar dimer model and the associ- 
ated random interface model. There has been a great deal of recent research 
on the dimer model but there are still many interesting open questions and 
new research avenues. In these lecture notes we will provide an introduction 
to the dimer model, leading up to results on limit shapes and fluctuations. 
There are a number of exercises and a few open questions at the end of each 
chapter. Our objective is not to give complete proofs, but we will at least 
attempt to indicate the ideas behind many of the proofs. Complete proofs 
can all be found in the papers in the bibliography. The main references we 
used for these notes are [TBI El El E] • The original papers of Kasteleyn [S] 
and Temperley/Fisher [27J are quite readable. 

Much previous work on the dimer model will not be discussed here; we'll 
focus here on dimers on planar, bipartite, periodic graphs. Of course one 
can study dimer models in which all or some of these assumptions have been 
weakened. While these more general models have been discussed and studied 
in the literature, there are fewer concrete results in this greater generality. It 
is our hope that these notes will inspire work on these more general problems. 

1.1 Dimer definitions 

A dimer covering, or perfect matching, of a graph is a subset of edges 
which covers every vertex exactly once, that is, every vertex is the endpoint 
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of exactly one edge. See Figure [TJ 

In these lectures we will deal only with bipartite planar graphs. A 
graph is bipartite when the vertices can be colored black and white in such 
a way that each edge connects vertices of different colors. Alternatively, 
this is equivalent to each cycle having even length. Kasteleyn showed how 
to enumerate the dimer covers of any planar graph, but the random sur- 
face interpretation we will discuss is valid only for bipartite graphs. There 
are many open problems involving dimer coverings of non-bipartite planar 
graphs, which at present we do not have tools to attack. However at present 
we have some nice tools to deal with periodic bipartite planar graphs. 

Our prototypical examples are the dimer models on Z 2 and the honey- 
comb graph. These are equivalent to, respectively, the domino tiling model 
(tilings with 2x1 rectangles) and the "lozenge tiling" model (tilings with 
60° rhombi) see Figures [1] and [2j 
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Figure 1: 



Dimers on the honeycomb have been studied in chemistry [TTI [7] where 
they are called Kekule structures. The honeycomb is after all the structure 
of graphite, each carbon atom sharing one double bond with a neighbor. 



1.2 Uniform random tilings 

Look at a larger domino picture and the lozenge picture, Figures [3] and 
SJ These are both uniform random tilings of the corresponding regions, 
that is, they are chosen from the distribution in which all tilings are equally 
weighted. In the first case there are about e 455 possible domino tilings and 
in the second, about e 1255 lozenge tilings Q These two pictures clearly display 
some very different behavior. The first picture appears homogeneous (and 

^^How do you pick a random sample from such a large space? 
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Figure 2: 



we'll prove that it is, in a well-defined sense), while in the second, the densities 
of the tiles of a given orientation vary throughout the region. A goal of these 
lectures is to understand this phenomenon, and indeed compute the limiting 
densities as well as other statistics of these and other tilings, in a setting of 
reasonably general boundary conditions. 

In figures and EJ we see a uniform random lozenge tiling of a triangular 
shape, and the same tiling rotated so that we see (with a little imagination) 
the fluctuations. These fluctuations are quite small, in fact of order ^/\o<gn 
for a similar triangle of side n. This picture should be compared with 
Figure [7] which shows the one-dimensional analog of the lozenge tiling of a 
triangle. It is just the graph of a simple random walk on Z of length n = 100 
conditioned to start and end at the origin. In this case the fluctuations 
are of order y/n. Indeed, if we rescale the vertical coordinate by \fn and 
the horizontal coordinate by n, the resulting curve converges as n — > oo to a 
Brownian bridge (a Brownian motion started at the origin and conditioned 
to return to the origin after time 1). 

The "scaling limit" of the fluctuations of the lozenge tiling of a triangle 
is a more complicated object, called the Gaussian free field. We can think of 
it as a Gaussian random function but in fact it is only a random distribution 
(weak function). We'll talk more about it later. 
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Figure 7: 



1.3 Limit shapes 

A lozenge tiling of a simply connected region is the projection along the 
(1, 1, Indirection of a piecewise linear surface in R 3 . This surface has pieces 
which are integer translates of the sides of the unit cube. Such a surface which 
projects injectively in the (1, 1, Indirection is called a stepped surface or 
skyscraper surface. Random lozenge tilings are therefore random stepped 
surfaces. A random lozenge tiling of a fixed region as in Figure H] is a random 
stepped surface spanning a fixed boundary curve in IR 3 . 

Domino tilings can also be interpreted as random tilings. Here the third 
coordinate is harder to visualize, but see Figure [8] for a definition. We'll 
see below that dimers on any bipartite graph can be interpreted as random 
surfaces. 

When thought of as random surfaces, these models have a limit shape 
phenomenon. This says that, for a fixed boundary curve in M 3 , or sequence 
of converging boundary curves in M 3 , if we take a random stepped surface on 
finer and finer lattices, then with probability tending to 1 the random surface 
will lie close to a fixed non-random surface, the so-called limit shape. So 
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Figure 8: The height is integer- valued on the faces and changes by ±1 across 
an unoccupied edge; the change is +1 when crossing so that a white vertex 
is on the left. 



this limit shape surface is not just the average surface but the only surface 
you will see if you take an extremely fine mesh... the measure is concentrating 
as the mesh size tends to zero on the delta-measure at this surface. 

Again the analogous property in one dimension is that the random curve 
of Figure [7] tends to a straight line when the mesh size tends to zero. (In 
order to see the fluctuations one needs to scale the two coordinates at dif- 
ferent rates. Here we are talking about scaling both coordinates by n: the 
fluctuations are zero on this scale.) 

As Okounkov puts it, the limit shape surface is in some sense the "most 
random" surface, in the sense that among all surfaces with the given bound- 
ary values, the limit shape is the one which has (overwhelmingly) the most 
discrete stepped surfaces nearby, in fact so many as to dwarf all other surfaces 
combined. 



1.4 Facets 

One thing to notice about lozenge tiling in figure H] is the presence of regions 
near the vertices of the hexagon where the lozenges are aligned. This phe- 
nomenon persists in the limit of small mesh size and in fact, in the limit 
shape surface there is a facet near each corner, where the limit shape is 
planar. This is a phenomenon which does not occur in one dimension (but 
see however [21]). The limit shapes for dimers generally contain facets and 
smooth (in fact analytic) curved regions separating these facets. In the facet 
the probability of a misaligned lozenge is zero in the limit of small mesh, 
and in fact one can show that the probability is exponentially small in the 
reciprocal of the mesh size. In particular the fluctuations away from the facet 



11 



are exponentially small. 
1.5 Measures 

What do we see if we zoom in to a point in figure HP That is, consider a 
sequence of such figures with the same fixed boundary but decreasing mesh 
size. Pick a point in the hexagon and consider the configuration restricted to 
a small window around that point, window which gets smaller as the mesh 
size goes to zero. One can imagine for example a window of side y/e when 
the mesh size is e. This gives a sequence of random tilings of (relative to the 
mesh size) larger and larger domains, and in the limit (assuming that a limit 
of these "local measures" exists) we will get a random tiling of the plane. 

We will see different types of behaviors, depending on which point we 
zoom in on. If we zoom in to a point in the facet, we will see in the limit a 
boring figure in which all tiles are aligned. This is an example of a measure on 
lozenge tilings of the plane which consists of a delta measure at a single tiling. 
This measure is an (uninteresting) example of an ergodic Gibbs measure (see 
definition below). If we zoom into a point in the non- frozen region, one can 
again ask what limiting measure on tilings of the plane is obtained. One 
of the important open problems is to understand this limiting measure, in 
particular to prove that the limit exists. Conjecturally it exists and only 
depends on the slope of the average surface at that point and not on any 
other property of the boundary conditions. For each possible slope (s, t) we 
will define below a measure fi Stt , and the local statistics conjecture states 
that, for any fixed boundary, fi s t is the measure which occurs in the limit 
at any point where the limit shape has slope (s,t). For certain boundary 
conditions this has been proved [T2] . 

These measures fi Sjt are discussed in the next section. 

1.5.1 Ergodic Gibbs Measures 

What are the natural probability measures on lozenge tilings of the whole 
plane? One might require for example that the measure be translation- 
invariant. In this case one can associate a slope or gradient to the measure, 
which is the expected change in height when you move in a given direction. 

Another natural condition to impose on a probability measure on tilings 
of the plane is that it be a Gibbs measure, that is, (in this case) a proba- 
bility measure which is the limit of the uniform measure on tilings of finite 
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regions, as the regions increase in size to fill out the whole plane. There is a 
generalization of this definition to dimers with edge weights. 

A remarkable theorem of Sheffield states that for each slope (s, t) for 
which there is an invariant measure, there is a unique ergodic Gibbs measure 
(ergodic means not a convex combination of other invariant measures). We 
denote this measure fj, Sjt . 

1.5.2 Phases 

Ergodic Gibbs measures come in three types, or phases, depending on the 
fluctuations of a typical (for that measure) surface. Suppose we fix the height 
at a face near the origin to be zero. A measure is said to be a frozen phase 
if the height fluctuations are finite almost surely, that is, the fluctuation of 
the surface away from its mean value is almost surely bounded, no matter 
how far away from the origin you are. A measure is said to be in a liquid 
phase if the fluctuations have variance increasing with increasing distance, 
that is, the variance of the height at a point tends to infinity almost surely 
for points farther and farther from the origin. Finally a measure is said to be 
a gaseous phase if the height fluctuations are unbounded but the variance 
of the height at a point is bounded independently of the distance from the 
origin. 

We'll see that for uniform lozenge or domino tilings we can have both 
liquid and frozen phases but not gaseous phases. An example of a graph 
for which we have a all three phases is the square octagon dimer model, 
see Figure O In general the classification of phases depends on algebraic 




Figure 9: 
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properties of the underlying graph and edge weights, as we'll see. 
1.6 Other random surface models 

There are a number of other statistical mechanical models with similar be- 
havior to the planar dimer model. The most well-known of these is the 
six-vertex model, or "square ice" model. This model consists of orientations 
of the edges of the square grid (or other planar graph of degree 4) with 
the restriction that at each vertex there are two "ingoing" and two "outgo- 
ing" arrows. This model also has an interpretation as a random interface 
model, with a limit shape phenomenon, facet formation, and so on. The 
dimer model is fundamentally easier to deal with than this model, and other 
models of this type, essentially because enumeration of dimer configurations 
can be accomplished using determinants, while the 6-vertex model requires 
more complicated and mysterious algebraic structures. Indeed, the quanti- 
ties which have been computed for these other models are been quite limited, 
and no one has succeeded in writing down a formula for a limit shape in any 
of these other models. 

2 The height function 

We show here how to define the height function for dimers on any bipartite 
graph. This allows us to give a "random surface" interpretation for dimers on 
any planar bipartite graph, using the height function as the third coordinate. 

2.1 Graph homology 

Let Q be a connected planar graph or connected graph embedded on a surface. 
Let A = A°(£?) be the space of functions on the vertices of Q. 

A flow, or 1-form, is a function uo on oriented edges of Q which is an- 
tisymmetric under changing orientation: if e is an edge connecting vertices 
V\ and v 2 and oriented from v\ to v 2 , then u>(e) = —u(—e) where by — e we 
denote the same edge oriented from v 2 to v±. Let A 1 be the space of 1-forms 
on Q. 

We define a linear operator d: A — > A 1 by dg(v\v 2 ) = g(v 2 ) — g{vi). The 
transpose of the operator d (for the natural bases indexed by vertices and 
edges) is d*, the divergence. The divergence d*u of a flow uj is a function 
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on vertices defined by 

d*cu(v) = u(e) 

e 

where the sum is over edges starting at e. A positive divergence at v means 
v is a source: water is flowing into the graph at v. A negative divergence is 
a sink: water is flowing out of the graph at v. 

We define a 2-form to be a function on oriented faces which is antisym- 
metric under changing orientation. Let A 2 be the space of 2-forms. We define 
d : A 1 — > A 2 , the curl, by: for a 1-form uo and face F, du(F) = ^ e uj(e), 
where the sum is over edges bounding F oriented consistently (cclw) with 
the orientation of F. The 1-forms cu for which dcu = are cocycles. The 
1-forms in the image of dA° are coboundaries, or gradient flows. 

Note that d o d as a map from A to A 2 is identically zero. Moreover for 
planar embeddings, the space of cocycles modulo the space of coboundaries 
is trivial, which is to say that every cocycle is a coboundary. Equivalently, 
the sequence 

A ° 4 A 1 4 A 2 - 

is exact. 

For a nonplanar graph embedded on a surface in such a way that every 
face is a topological disk, the space of cocycles modulo the space of cobound- 
aries is the 1-homology of the surface. In particular for a graph embedded in 
this way on a torus, the 1-homology is M 2 . The 1-forms which are nontrivial 
in homology are those which have non-zero net flux across the horizontal 
and/or vertical loops around the torus. The point in M 2 whose coordinates 
give these two net fluxes is the homology class of the 1-form. 

2.2 Heights 

Given a dimer cover M of a bipartite graph, there is a naturally associated 
flow ui 'm'- it has value 1 on each edge occupied by a dimer, when the edge is 
oriented from its white vertex to its black vertex. Other edges have flow 0. 
This flow has divergence d*UM = ±1 at white, respectively black vertices. 

In particular given two dimer covers of the same graph, the difference of 
their flows is a divergence-free flow. 

Divergence-free flows on planar graphs are dual to gradient flows, that is, 
Given a divergence- free flow and a fixed face /o, one can define a function h 
on all faces of Q as follows: h(fo) = 0, and for any other face /, h(f) is the 
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net flow crossing (from left to right) a path in the dual graph from f to /. 
The fact that the flow is divergence-free means that h(f) does not depend 
on the path from / to /. 

To define the height function for a dimer cover, we proceed as follows. 
Fix a flow ujq with divergence 1 at white vertices and divergence —1 at black 
vertices (ujq might come from a fixed dimer cover but this is not necessary). 
Now for any dimer cover M, let % be its corresponding flow. Then the dif- 
ference Um — Wq is a divergence-free flow. Let h = h(M) be the corresponding 
function on faces of Q (starting from some fixed face /o). Then h(M) is the 
height function of M. It is a function on faces of Q. Note that if Mi, M2 
are two dimer coverings then h{M\) — h(M 2 ) does not depend on the choice 
of u . In particular ujq is just a "base point" for h and the more natural 
quantity is really the difference h(Mi) — h{M2). 

For lozenge tilings a natural choice of base point flow is the flow ujq with 
value 1/3 on each edge (oriented from white to black). In this case it is not 
hard to see that h(M) is just, up to an additive constant, the height (distance 
from the plane x + y + z = 0) when we think of a lozenge tiling as a stepped 
surface. Another base flow, u a , which is one we make the most use of below, 
is the flow which is +1 on all horizontal edges and zero on other edges. For 
this base flow the three axis planes (that is, the dimer covers using all edges 
of one orientation) have slopes (0, 0), (1, 0) and (0, 1). 

Note that for a bounded subgraph in the honeycomb, the flow u or uo a will 
typically not have divergence ±1 at white/black vertices on the boundary. 
This just means that the boundary height function is nonzero: thought of as 
a lozenge tiling, the boundary curve is not flat. 

Similarly, for domino tilings a natural choice for u$ is 1/4 on each edge, 
oriented from white to black. This is the definition of the height function 
for dominos as defined by Thurston [2H] and also (1/4 of) that in Figure El 
Thurston used it to give a linear-time algorithm for tiling simply-connected 
planar regions with dominos (that is deciding whether or not a tiling exists 
and building one if there is one). 

For a nonplanar graph embedded on a surface, the flows can be defined 
as above, but not the height functions in general. This is because there 
might be some period, or height change, along topologically nontrivial closed 
dual loops. This period is however a homological invariant in the sense 
that two homologous loops will have the same period. For example for a 
graph embedded on a torus, we can define two periods h x , h y for loops in 
homology classes (1, 0) and (0, 1) respectively. Then any loop in homology 
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class (x, y) G Z 2 will have height change xh x + yh y . 

Exercise 1. From an 8 x 8 checkerboard, two squares are removed from 
opposite corners. In how many ways can the resulting figure be tiled with 
dominos? 

Exercise 2. Using the height function defined in Figure 0, find the lowest 
and highest tilings of a 2m x 2n rectangle. 

Exercise 3. Take two hexagons with edge lengths 11, 13, 11, 13, 11, 13 in 
cyclic order, and glue them together along one of their sides of length 11. 
Can the resulting figure be tiled with lozenges? What if we glue along a side 
of length 13? 

Exercise 4. From a 2n x 2n square, a subset of the squares on the lower 
edge are removed. For which subsets does there exist a domino tiling of the 
resulting region? 

Exercise 5. For a finite graph Q, let Q be the set of flows uj satifying u>(e) e 
[0, 1] when e is oriented from white to black, and having divergence ±1 and 
white/black vertices, respectively. Show that Q is a convex polytope. Show 
that the dimer coverings of Q are the vertices of this polytope. In particular 
ifQ^tf) then there is a dimer cover. 

3 Kasteleyn theory 

We show how to compute the number of dimer coverings of any bipartite 
planar graph using the KTF (Kasteleyn- Temperley-Fisher) technique. While 
this technique extends to nonbipartite planar graphs, we will have no use for 
this generality here. 

3.1 The Boltzmann measure 

To a finite bipartite planar graph Q = (BUW,E) with a positive real function 
w : E — > 1R >0 on edges, we define a probablity measure \i = fi(Q, w) on dimer 
covers by, for a dimer covering M, 

MM) = i J] w(e) 

eeM 
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where the sum is over all dimers of M and where Z is a normalization constant 
called the partition function, defined to be 

Z = }ZU 

M e<=M 

3.2 Gauge equivalence 

If we change the weight function w by multiplying edge weights of all edges 
incident to a single vertex v by a constant A, then Z is multiplied by A, and 
the measure \x does not change, since each dimer cover uses exactly one of 
these edges. So we define two weight functions w, w' to be equivalent, w ~ w', 
if one can be obtained from the other by a sequence of such multiplications. 
It is not hard to show that w ~ w' if and only if the alternating products 
along faces are equal: given a face with edges e±, e 2 , . . . , e 2k in cyclic order, 
the quotients 

w(ei)w(e 3 ) . . . w(e 2k -i) 
w(e 2 )w(e 4 ) . ..w(e 2k ) 

and 

w'(ei)w'(e 3 ) . . . w'(e 2k -i) 
w'(e 2 )w'(e A ) . ..w'(e 2k ) 

(which we call alternating products) must be equal. For nonplanar graphs, 
gauge equivalence is the same as having equal alternating products along all 
cycles (and it suffices to consider cycles in a homology basis). 

The proof is easy, and in fact is just a little homology theory: logw is 
a 1-form (if you orient edges from white to black) and w is equivalent to 
w' if and only if logw — log«/ = df for some function / G A . However 
logw — logw' = df for a planar graph is equivalent to dlogw = dlogw' 
which can be interpreted as saying that the alternating products of w and 
w' on cycles are equal. 

3.3 Kasteleyn weighting 

A Kasteleyn weighting of a planar bipartite graph is a choice of sign for 
each undirected edge with the property that each face with mod 4 edges 
has an odd number of — signs and each face with 2 mod 4 edges has an even 
number of — signs. 
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In certain circumstances it will be convenient to use complex numbers of 
modulus 1 rather than signs ±1. In this case the condition is that the alter- 
nating product of edge weights (as defined above) around a face is negative 
real or positive real depending on whether the face has or 2 mod 4 edges. 

This condition appears mysterious at first but we'll see why it is impor- 
tant below^. We can see as in the previous section that any two Kasteleyn 
weightings are gauge equivalent: they can be obtained one from the other by 
a sequence of operations consisting of multiplying all edges at a vertex by a 
constant. 

The existence of a Kasteleyn weighting is also easily established for ex- 
ample using spanning trees. We leave this fact to the reader, as well as the 
proof of the following (easily proved by induction) 

Lemma 1. Given a cycle of length 2k enclosing I points, the alternating 
product of signs around this cycle is (— l) 1+fc+£ . 

Note finally that for the (edge-weighted) honeycomb graph, all faces have 
2 mod 4 edges and so no signs are necessary in the Kasteleyn weighting. 

3.4 Kasteleyn matrix 

A Kasteleyn matrix is a weighted, signed adjacency matrix of the graph Q. 
Given a Kasteleyn weighting of Q, define a \B \ x \W\ matrix K by K(b, w) = 
if there is no edge from w to b, otherwise K(b, w) is the Kasteleyn weighting 
times the edge weight u'(bw). 

For the graph in Figure [10] with Kasteleyn weighting indicated, the Kaste- 
leyn matrix is 

a 1 
1 -b 1 
1 c 

Note that gauge transformation corresponds to pre- or post-multiplication 
of if by a diagonal matrix. 

Theorem 2 (0E7]). Z = | detK\. 

In the example, the determinant is —a — c — abc. 

2 The condition might appear more natural if we note that the alternating product is 
required to be e^ lN ^ 2 where N is the number of triangles in a triangulation of the face 
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Proof. If K is not square the determinant is zero and there are no dimer 
coverings (each dimer covers one white and one black vertex). If K is a 
square n x n matrix, We expand 

detK — ^ s g n (°")^( b i> w CT(1 )) . . . K(b n , w^ij). (1) 

cr£S n 

Each term is zero unless it pairs each black vertex with a unique neighboring 
white vertex. So there is one term for each dimer covering, and the modulus 
of this term is the product of its edge weights. We need only check that the 
signs of the nonzero terms are all equal. 

Let us compare the signs of two different nonzero terms. Given two 
dimer coverings, we can draw them simultaneously on Q. We get a set of 
doubled edges and loops. To convert one dimer covering into the other, we 
can take a loop and move every second dimer (that is, dimer from the first 
covering) cyclically around by one edge so that they match the dimers from 
the second covering. When we do this operation for a single loop of length 
2k, we are changing the permutation a by a fc-cycle. Note that by Lemma [1] 
above the sign change of the edge weights in the corresponding term in (jTJ 
is ±1 depending on whether 2k is 2 mod 4 or mod 4 (since i is even there), 
exactly the same sign change as occurs in sgn(tx). These two sign changes 
cancel, showing that these two coverings (and hence any two coverings) have 
the same sign. □ 

A simpler proof of this theorem for honeycomb graphs — which avoids 
using Lemma [1] — goes as follows: if two dimer coverings differ only on a 
single face, that is, an operation of the type in Figure [TT] converts one cover 
into the other, then these coverings have the same sign in the expansion of 
the determinant, because the hexagon flip changes a by a 3-cycle which is 
an even permutation. Thus it suffices to notice that any two dimer coverings 
can be obtained from one another by a sequence of hexagon flips. This 
can be seen using the lozenge tiling picture since applying a hexagon flip is 
equivalent to adding or subtracting a cube from the stepped surface. Any 
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Figure 11: 



two surfaces with the same connected boundary can be obtained from one 
another by adding and/or subtracting cubes. 

While there is a version of Theorem [2] (using Pfaffians) for non-bipartite 
planar graphs, there is no corresponding sign trick for nonplanar graphs in 
general (the exact condition is that a graph has a Kasteleyn weighting if and 
only if it does not have -fQ^ as minor [18J). 

3.5 Local statistics 

There is an important corollary to Theorem [2J 

Corollary 3 ([TT]). Given a set of edges X = {wibi, . . . , w^b^}, the proba- 
bility that all edges in X occur in a dimer cover is 

k \ 

JJi^b^Wj) J det(K 1 (w i ,b j ))i< itj < k . 
i=i / 

The proof uses the Jacobi Lemma that says that a minor of a matrix A 
is detA times the complementary minor of A" 1 . 

The advantage of this result is that the probability of a set of k edges 
being present is only a k x k determinant, independently of the size of the 
graph. One needs only be able to compute K~ l . In fact the corollary is valid 
even for infinite graphs, once K~ x has been appropriately defined. 

Corollary [3] shows that the edges form a determinantal (point) pro- 
cess. Such a process is defined by the fact that the probability of a set 
of k "points" pi, . . . ,pk is a determinant of a k x k matrix M with entries 
M(pi,pj). Here the points in the process are edges of Q and M(pi,pj) = 
K(bi, Wj)i^ _1 (wj, bj) where Wj is the white vertex of edge Pi and bj is the 
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black vertex of edge pj. See [22] for an introduction to determinantal pro- 
cesses. 

Exercise 6. Without using Kasteleyn theory, compute the number of domino 
tilings of a 2 x n and 3 x n rectangle. 

Exercise 7. A classical combinatorial result of Macmahon says that the num- 
ber of lozenge tilings of an A x B x C hexagon ( that is, hexagon with edges 
A, B, C, A, B, C in cyclic order) is 

A B C . + 

nnri i+J+fc _2 

i=i j=i k=i j 

What is the probability that, in a uniform random tiling, there are two lozenges 
adjacent to a chosen corner of an n x n x n hexagon^? 

Exercise 8. Show that a planar bipartite graph has a Kasteleyn weighting 
and prove LemmaUl 

Exercise 9. Verify Corollary^ for the example in Figure UR 



4 Partition function 

4.1 Rectangle 

Here is the simplest example. Assume mn is even and let Q m ,n be the m x n 
square grid. Its vertices are V = {1,2,..., m} x {1,2, ... ,n} and edges 
connect nearest neighbors. Let Z m ^ n be the partition function for dimers 
with edge weights 1. This is just the number of dimer coverings of Q m ,n- 

A Kasteleyn weighting is obtained by putting weight 1 on horizontal edges 
and i = a/— T on vertical edges. Since each face has four edges, the condition 
in (13.31) is satisfied^ 

3 In fact a more general formula holds: If we weight configurations with g vok,mc j then 

A B c I _ qi+j+k-l 

z = nn n l _ ^+^-2 ■ 

4 A weighting gauge equivalent to this one, and using only weights ±1, is to weight 
alternate columns of vertical edges by —1 and all other edge +1. This was the weighting 
originally used by Kasteleyn [Sj; our current weighting (introduced by Percus [53]) is 
slightly easier for our purposes. 



22 



The corresponding Kasteleyn matrix K is an mn/2 x mn/2 matrix (recall 

f K 

that K is a | W | x \B\ matrix). The eigenvalues of the matrix K 



K l 

are in fact simpler to compute. Let z = e l7V 3/( rn + 1 ) and w = e ; mk /i n + l ) , Then 
the function 

f jth (x, y) = (z* - z-*){ W y - w-y) = -4 sin(^-) rinA 

is an eigenvector of K with eigenvalue z + ^ + i (w + —). To see this, check 
that 

A/(ar, 2/) = /(s + 1,?/) + /(a; - 1, y) + y + 1) + if(x, y-l) 

when (x,y) is not on the boundary of Q, and also true when / is on the 
boundary assuming we extend / to be zero just outside the boundary, i.e. 
when x — Oovy — Oovx — m + lovy — n+1. 

As j, k vary in (l,m) x (l,n) the eigenfunctions fj^ are independent (a 
well-known fact from Fourier series). Therefore we have a complete diago- 
nalization of the matrix K, leading to 

rrrf2cos^- + 2zcos^- . (2) 

m+l n+1 

j=ik=i J 

Here the square root comes from the fact that det K = (det K) 2 . 

Note that if m, n are both odd then this expression is zero because of the 
term j = (m + 1) /2 and k = (n + l)/2 in (J2J. 

For example Z 8)8 = 12988816. For large m, n we have 

1 y r r 

lim log Z m n = — — / / log(2 cos 6 + 2i cos 

m,n^oo mn 2rc z J J 



which can be shown to be equal to G/n, where G is Catalan's constant 
ri_i_i i l_ 

— 1 3 2 f 52 

Exercise 10. STiotu i/iis. 



23 



4.2 Torus 



A graph on a torus does not in general have a Kasteleyn weighting. However 
we can still make a "local" Kasteleyn matrix whose determinant can be used 
to count dimer covers. 

Rather than show this in general, let us work out a detailed example. Let 
H n be the honeycomb lattice on a torus, as in Figure [121 which shows H$. 
It has n 2 black vertices and n 2 white vertices, and 3n 2 edges. Weight edges 




Figure 12: Fundamental domain for the honeycomb graph on a torus. 

a, b, c according to direction: a the horizontal, b the NW-SE edges and c the 
NE-SW edges. Let x and y be the directions indicated, and h x and h y be the 
height changes along the paths winding around the torus in these directions, 
defined using as base flow uo a (the flow coming from the "all-a" dimer cover). 
In particular, h x , h y are the number of b-, respectively c-type edges crossed 
on a path winding around the torus in the x, respectively y direction. 
We have the following lemma. 

Lemma 4. The height change (h x ,h y ) of a dimer covering of H n is deter- 
mined by the numbers N a , Nf,, N c of dimers of each type a, b, c, as follows: 
h x = ^ and h v = 

Proof. We can compute the height change h x along any path winding around 
the torus in the x-direction. Take for example the dashed path in Figure [12j 
The height change h x is the number of dimers crossing this path. We get the 
same value h x for any of the n translates of this path to other y coordinates. 
Summing over all translates gives nh x = Nf, and dividing by n gives the 
result. The same argument holds in the y-direction. □ 
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Let K n be the weighted adjacency matrix of H n . That is K(b,w) = if 
there is no edge from b to w and otherwise K(b, w) = a, b or c according to 
direction. 

From the proof of Theorem [2] we can see that det K is a weighted, signed 
sum of dimer coverings. Our next goal is to determine the signs. 

Lemma 5. The sign of a dimer covering in det K depends only on its height 
change (h x , h y ) modulo 2. Three of the four parity classes gives the same sign 
and the fourth has the opposite sign. 

Proof. Let iV 6 , N c be the number of b and c type edges in a cover. If we take 
the union of a covering with the covering consisting of all a-type edges, we 
can compute the sign of the cover by the product of the sign changes when 
shifting along each loop. The number of loops is q = GCF(h x , h y ) and each of 
these has homology class (h y /q, h x /q) (note that the b edges contribute to h x 
but to the (0, 1) homology class). The length of each loop is ^(h x + h y ) and 
so each loop contributes sign (-i) 1+ ^ h * +h y\ for a total of (-l^+K^+M. 
Note that q is even if and only if h x and h y are both even. So if n is even 
then the sign is —1 unless (h x , h y ) = (0, 0) mod 2. If n is odd the sign is +1 
unless (h x , h y ) = (1, 1) mod 2. □ 

In particular det K, when expanded as a polynomial in a, b and c, has co- 
efficients which count coverings with particular height changes. For example 
for n odd we can write 

detK = J2(-l) hxhv C hx)hy a n(n ~ hx ~ hv) b nhx c nh y 

since h x h y is odd exactly when h x , h y are both odd. 
Define Z Q0 = Z 00 (a, b, c) to be this expression. Define 

Z 10 (a,b,c) =Z QQ (a,be iri/n ,c), 

Z 01 (a,b,c)=Z 00 (a,b,ce m / n ), 

and 

Z 11 (a,b,c) = Z 00 (a,be^ n ,ce^ n ). 
Then one can verify that, when n is odd, 

Z = -(Z 00 + Ziq + Z i — Zu) (3) 

which is equivalent to Kasteleyn's expression [8] for the partition function. 
The case n even is similar and left to the reader. 



25 



4.3 Partition function 



Dealing with a torus makes computing the determinant much easier, since the 
graph now has many translational symmetries. The matrix K commutes with 
translation operators and so can be simultaneously diagonalized with them. 
(In other words, we can use a Fourier basis.) Simultaneous eigenfunctions of 
the horizontal and vertical translations (and K) are the exponential functions 
fz,w{%, y) — z~ x w~ y where z n = 1 = w n . There are n choices for z and n for 
w leading to a complete set of eigenfunctions. The corresponding eigenvalue 
for K is a + bz + cw. 
In particular 



det K = Z 00 = Yl ] [ a + bz + 



cw. 

z n—i w n—i 



This leads to 



Z w = JJ JJ a + bz + , 

z n =-l w n =l 

Z i=Y[ n a + bz + , 

z n =l w n =—l 

Z n = |[ ]"[ a + bz + 



cw. 

z n =-l w n =-l 



4.4 Height change distribution 

From the expression for Z 00 one can try to estimate, for a,b,c — 1, the size 
of the various coefficients. Assuming that n is a multiple of 3, the largest 
coefficient Ch x ,h y occurs when (h x ,h y ) = (n/3, n/3), that is, (N a , ]V&, N c ) = 
(T' T' t)' Boutillier and de Tiliere [2j showed that, letting C = C n / 3jn /3 
and (/ix,/^) = (j + f,/c + f), 

C^=^e- c ^' fc+fc2 )(l + 0(l/n)) 

where c is a constant. So for large n the height change of a uniform ran- 
dom tiling on a torus has a distribution converging to a discrete Gaussian 
distribution centered at (n/3, n/3). 

Exercise 11. VF/iai szgns should we put in (Tj|) in i/ie case n is even? 

Exercise 12. Compute the partition function for dinners on a square grid on 
a cylinder of circumference 4n + 2 and height m. Weight horizontal edges a 
and vertical edges b. 
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5 Gibbs measures 



5.1 Definition 

Let X = X(Q) be the set of dimer coverings of a graph Q, possibly infinite, 
with edge weight function w. Recall the definition of the Boltzmann prob- 
ability measure on X{Q) when Q is finite: a dimer covering has probability 
proportional to the product of its edge weights. When Q is infinite, this defi- 
nition will of course not work. For an infinite graph Q, a probability measure 
on X is a Gibbs measure if it is a weak limit of Boltzmann measures on a 
sequence of finite subgraphs of Q filling out Q. By this we mean, for any finite 
subset of Q, the probability of any particular configuration occurring on this 
finite set converges. That is, the probabilities of cylinder sets converge. 
Here a cylinder set is a subset of X(Q) consisting of all coverings which have 
a particular configuration on a fixed set of vertices. For example the set of 
coverings containing a particular edge is a cylinder set. 

For a sequence of Boltzmann measures on increasing graphs, the limiting 
measure may not exist, but subsequential limits will always exist. The limit 
may not be unique, however; that is, it may depend on the approximating 
sequence of finite graphs. This will be the case for dimers and it is this 
non-uniqueness which makes the dimer model interesting. 

The important property of Gibbs measures is the following. Let A be a 
cylinder set defined by the presence of a finite set of edges. Let B be another 
cylinder set defined by a different set of edges but using the same set of 
vertices. Then, for the approximating Boltzmann measures, the ratio of the 
measures of these two cylinder sets is equal to the ratio of the product of the 
edge weights in A and the product of the edge weights in B. This is true along 
the finite growing sequence of graphs and so in particular the same is true 
for the limiting Gibbs measures. In fact this property characterizes Gibbs 
measures: given a finite set of vertices, the measure on this set of vertices 
conditioned on the exterior (that is, integrating over the configuration in the 
exterior) is just the Boltzmann measure on this finite set. 

5.2 Periodic graphs 

We will be interested in the case when Q is a periodic bipartite planar 

graph, that is, a planar bipartite weighted graph on which translations in Z 2 
(or some other rank-2 lattice) act by weight-preserving and color-preserving 
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isomorphisms. Here by color-preserving isomorphisms, we mean, isomor- 
phism which maps white vertices to white and black vertices to black. Note 
for example that for the graph Q = Z 2 with nearest neighbor edges, the lat- 
tice generated by (2, 0) and (1, 1) acts by color-preserving isomorphisms, but 
Z 2 itself does not act by color-preserving isomorphisms. So the fundamental 
domain contains two vertices, one white and one black. 

For simplicity we will assume our periodic graphs are embedded so that 
the lattice of weight- and color-preserving isomorphisms is Z 2 , so that we can 
describe a translation using a pair of integers. 



5.3 Ergodic Gibbs measures 

For a periodic graph Q, a translation-invariant measure on X(Q) is sim- 
ply one for which the measure of a subset of X{Q) in invariant under the 
translation-isomorphism action. 

The slope (s, t) of a translation-invariant measure is the expected height 
change in the (1,0) and (0,1) directions, that is, s is the expected height 
change between a face and its translate by (1,0) and t is the expected height 
change between a face and its translate by (0, 1). 

A ergodic Gibbs measure, or EGM, is one in which translation-invariant 
sets have measure or 1. Typical examples of translation-invariant sets are: 
the set of coverings which contain a translate of a particular pattern. 

Theorem 6 (Sheffield [25J). For the dimer model on a periodic planar bipar- 
tite periodically edge-weighted graph, for each slope (s, t) for which there exists 
a translation-invariant measure, there exists a unique EGM fi Sjt - Moreover 
every EGM is of this type for some s, t. 

In particular we can classify EGMs by their slopes. 

The existence is not hard to establish by taking limits of Boltzmann 
measures on larger and larger tori with restricted height changes (h x ,h y ), see 
below. The uniqueness is much harder; we won't discuss this here. 



5.4 Constructing EGMs 

Going back to our torus H n , note that the number of a, b, and c type edges 
are multiples of n and satisfy 

N a N b N 
< — , — , — < n 

n n n 
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N a N b t N c 

1 1 —n 

n n n 

and in fact it is not hard to see that every triple of integers in this triangle 
can occur. Recalling that the height changes were related to these by h x = 
Nb/n, h y = N c /n, and the average slope (s,t) is defined by (h x /n, h y /n), we 
have that (s, t) lies in the triangle 

{(s,t) : < s,t,s + t < 1}. 

We conclude from Theorem [6] that there is a unique EGM for every (s,t) 
in this triangle. We denote this EGM /i s t . 

We can construct fi S}t as follows. On X(H n ) restrict to configurations for 
which (h x , h y ) = ([nsj, |_ n ^J)- Let jj^} be the Boltzmann measure on X(H n ) 
conditioned on this set. Any limit of the /4t as n — > oo is a translation- 
invariant Gibbs measure of slope (s,t). Ergodicity follows from uniqueness: 
the set of translation-invariant Gibbs measures of slope (s, t) is convex and 
its extreme points are the ergodic ones; since there is a unique ergodic one 
this convex set must be reduced to a point. 

In practice, it is hard to apply this construction, since conditioning is a 
tricky business. The next section gives an alternate construction. 

5.5 Magnetic field 

The weights a, b, c are playing dual roles. On the one hand they are variables 
in the partition function whose exponents determine the height change of a 
covering. On the other hand by putting positive real values in for a, b, c we 
reweight the different coverings. This reweighting has the property that it 
depends only on h x ,h y , that is, two configurations with the same h x , h y are 
reweighted by the same quantity. As a consequence putting in weights a, b, c 
has the effect of changing the average value of (h x , h y ) for a random dimer 
cover of T n . However a random dimer cover of a planar region is unaffected 
by this reweighting since the Boltzmann measures are unchanged. 

There is in fact a law of large numbers for covers of T n : as the torus gets 
large, the slope of a random tiling (with edges weights a, b, c) is concentrating 
on a fixed value (s, t), where s, t are functions of a, b, c. We'll see this below. 

This reweighting is analogous to performing a simple random walk in 1 
dimension using a biased coin. The drift of the random walk is a function of 
the bias of the coin. In our case we can think of a, b, c as a bias which affects 
the average slope (which corresponds to the drift). 
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We computed Z n (a, b, c) above; one can find (s, t) as a function of (a, b, c) 
from this formula. It is easier to use the asymptotic expression for Z which 
we compute in the next section. 

Exercise 13. For (s, t) on the boundary of the triangle of possible slopes, 
describe the corresponding measure // s ^. (Hint: try the corners first.) 



6 Uniform honeycomb dimers 

Recall the expressions 

Z 00 = ]~| Y\ a + bz + cw, 



z n =l w n =l 



z n = -l w n =l 

Z ol = ] [ [ f a + bz + cw, 

z n =l w n =—l 

Z n = Y\ JJ a + bz + cw, 



-1 w n =-l 



and (for n odd) 



Z — -{Zqo + Zio + Z i — Z\\). 



When n is large, we can estimate these quantities using integrals. Indeed, 
the logs of the right-hand sides are different Riemann integrals for log (a + 
bz + cw). Note that 



max Z ar < Z < 2 max Z ar , 

<7,T a.T 



so the nth root of Z and the nth root of the maximum of the Z UT have the 
same limits. 

It requires a bit of work [5] to show that 



Theorem 7. 



1- 1 , ry 1 f f T / 7 .dzdw 

hm — logZ = — — — ■ / / Logfa + bz + cw) . 

™) is 1 Js 1 z w 



1 

n^oo n 2 to (27T2 
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Here Log denotes the principal branch. By symmetry under complex 
conjugation the integral is real. This normalized logarithm of the partition 
function is called the free energy F. 

The difficulty in this theorem is that the integrand has two singularities 
(assuming a, b, c satisfy the triangle inequality) and the Riemann sums may 
become very small if a point falls near one of these singularities^]. However 
since the four integrals are Riemann sums on staggered lattices, at most one 
of the four can become small, (and it will be come small at both singularities 
by symmetry under complex conjugation.) Therefore either three or all four 
of the Riemann sums are actually good approximations to the integral. This 
is enough to show that log Z, when normalized, converges to the integral. 

This type of integral is well-studied. It is known at the Mahler measure 
of the polynomial P(z, w) = a + bz + cw. We'll evaluate it for general a, b, c 
below. Note that when a > b + c the integral can be evaluated quickly by 
residues and gives F = log a. Similarly when b > a + c or c > a + b we get 
F = log b or F = log c respectively. 



6.1 Inverse Kasteleyn matrix 

To compute probabilities of certain edges occurring in a dimer covering of 
a torus, one needs a linear combination of minors of four inverse Kasteleyn 
matrices. Again there are some difficulties due to the presence of zeros in 
the integrand, see [To] . In the limit n — > oo, however, there is a simple 
expression for the multiple edge probabilities. It is identical to the statement 
in Corollary [31 except that we must use the infinite matrix K~ l defined by: 

, , . 1 f z~ v w x dz dw 

K (W ,0,b x , y ) = — — r-r / . 

(2m) z J a + bz + cw z w 

Here we are using special coordinates for the vertices; wo,o is the white vertex 
at the origin and h x>y corresponds to a black vertex at ei +x(e 3 — ei) +y(e\ — 
e 2 ), where ei,e 2 ,e 3 are the unit vectors in the directions of the three cube 
roots of 1. 

As an example, suppose that a, b, c satisfy the triangle inequality. Let 
a , Ob, C be the angles opposite sides a, b, c in a Euclidean triangle. The 

5 To see where the singularities occur, make a triangle with edge lengths a, b, c; think 
about it sitting in the complex plane with edges a,bz,cw where \z\ — \w\ = 1. The two 
possible orientations of triangle give complex conjugate solutions (z, w), (z,w). 
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probability of the a edge w 0j obo,o i n a dimer covering of the honeycomb is 



Pr(w , b ,o) = aK 1 (w ,o,b 0) o) = 

Doing a contour integral over w gives 

1 



1 



47T 5 



S 1 JS 



dz dw 



i a + bz + cw iz iw 



which simplifies to 



2m 



1 



a dz 

a+bz\>c a ~\~bz Z 



bz 



,i(-rr-e c ) 



-i(n-Bc) 



a + bz 

= 0_a 

7T 

Note that this tends to 1 or when the triangle degenerates, that is, when 
one of a, b, c exceeds the sum of the other two. This indicates that the dimer 
covering becomes frozen: when a > b + c only a- type edges are present. 



6.2 Decay of correlations 

One might suspect that dimers which are far apart are uncorrelated: that 
is, the joint probability in a random dimer cover of two (or more) edges 
which are far apart, is close to the product of their individual probabili- 
ties. This is indeed the case (unless s, t is on the boundary of the triangle 
of allowed slopes), and the error, or correlation, defined by Pr(ei&e2) — 
Pr(ei) Pr(e2), is an important measure of how quickly information is lost 
with distance in the covering. By Corollary [31 this correlation is a constant 
times K~ 1 (bi,w 2 )K~ 1 (b 2 ,wi) if the two edges are biWx and b 2 w 2 . 

The values of i^ _1 (w 0i o, b^) are Fourier coefficients of (a + bz + cw)^ 1 . 
When b and w are far from each other, what can we say about fC _1 (w, b)? 
The Fourier coefficients of an analytic function on the torus {\z\ = \w\ = 1} 
decay exponentially fast. However when a, b, c satisfy the triangle inequality, 
the function (a + bz + cw)' 1 is not analytic; it has two simple poles (as we 
discussed earlier) on {\z\ = \w\ = 1}. 

The size of the Fourier coefficients is governed by the behavior at these 
poles, in fact the Fourier coefficients decay linearly in \x\ + \y\. As a conse- 
quence the correlation of distant edges decays quadratically in the distance 
between them. 
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This is an important observation which we use later as well: the polyno- 
mial decay of correlations in the dimer model is a consequence of the zeros 
of a + bz + cw on the torus. In more general situations we'll have a different 
polynomial P(z,w), and it will be important to find out where its zeros lie 
on the unit torus. 

6.3 Height fluctuations 

As in section 16.11 above we can similarly compute, for k ^ 0, 

—if u \ - sin(k6 b ) 



K (w , ,bi„ 



ixkb 



By symmetry (or maybe there's an easy way to evaluate this directly, I don't 
know) for k ^ 0, 

K (w ,o,b fcffe ) = . 

Let a k = _ sin (^ 6> ") for A; ^ and a = 9 a /n. 

Given any set of k edges of type a in the vertical column passing through 
the edge w ,obo,o, say the edges w n ^ nj .b n ^ for j = 1, ... k, by Corollary [3] the 
probability of these k edges all occurring simultaneously is det(a n ^_ nj )i<j ) -,<fc. 

So the presence of the edges in this column forms a determinantal process 
with kernel My = a^j. Determinantal processes [26J have the property that 
the number of points in any region is a sum of independent (but not necessar- 
ily identically distributed) Bernoulli random variables. The random variables 
have probabilities which are the eigenvalues of the kernel M restricted to the 
domain in question. In particular if the variance in the number of points in 
an interval is large, this number is approximately Gaussian. 

We can compute the variance in the number of points in an interval of 
length k as follows. Let be the k x k matrix M^j = a^j for 1 < i, j < k. 
The sum of the variances of the Bernoullis is just the sum of A(l — A) over 
the eigenvalues of This is just Tr(Mfc(Jfc — M^)) (where is the identity 
matrix). With a bit of work (and the Fourier transform of the function 
f(9) = \9\ (n — \9\) for 9 G [— ir, 7r]) one arrives at 

Tr(M k (I - M k )) = ka {l - a ) - aj(2k - 2) - a\(2k - 4) a\_ x (4) 

= — logfc + 0(l). 

7l z 
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We conclude that the variance in the height difference between the face 
at the origin and the face at (k, k), that is, k lattice spacings vertically away 
from the origin, is proportional to the log of the distance between the faces. 
In particular this allows us to conclude that the height difference between 
these points tends to a Gaussian when the points are far apart. 

A similar argument gives the same height difference distribution for any 
two faces at distance k (up to lower order terms). 

Exercise 14. For the uniform measure a = b = c = 1, compute the proba- 
bility that the face at the origin has three of its edges matched, that is, looks 
like one of the two configurations in Figure [771 

Exercise 15. Compute the constant term in the expression Q). 

Exercise 16. How would you modify to compute the expected parity of 
the height change from the face at the origin to the face at (k, k) ? 

7 Legendre duality 

Above we computed the partition function for dimer coverings of the torus 
with edge weights a, b, c. Here we would like to compute the number of 
dimer coverings of a torus with uniform weights but with fixed slope (s, t). 
Surprisingly, these computations are closely related. Indeed, we saw that, 
with edge weights a, b, c, all coverings with height change h x , h y have the 
same weight a nhz b nhx c nhy (we defined h z to be h z = n — h x — h y ). 

So we just need to extract the coefficient Ch x ,h y in the expansion of the 
partition function 



This can be done as follows. First, choose positive reals a, b, c (if we can) 
so that the term Ch x ,h y < inhz b nh:c c nhv is the largest term. If we're lucky, this 
term and terms with nearby h x , h y dominate the sum, in the sense that all 
the other terms add up to a negligible amount compared to these terms. In 
that case we can use the estimate 



Z n (a,b,c)= Y, C h*,h v °< nh *b nh *c nh *. 
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and so 

C hxAy « Z n (a, b } c)a- n ^ l - s -%- n2s c- n2t . 

\h x /n—s\<€,\h y /n—t\<€ 

These estimates can in fact be made rigorous. One needs only check that 
for fixed edge weights a,b,c, as n — > oo the height change concentrates on 
a fixed value (s,t), that is, the variance in the average slope (s,t) tends to 
zero. 

We can conclude that the growth rate (which we denote —a(s,t)) of the 
number of stepped surfaces of fixed slope (s, t) is 



- a(s, t) = lim \ log C, 



n— >oo fl 



h x ,h v 



lim — (\ogZ n — n 2 (l — s — t) log a — n 2 t\ogb — n 2 slogc) 



= \ogZ - p a \oga - p b \ogb - p c \ogc, (5) 

where (p a ,Pb,Pc) = (1 — s — t,s,t) are the probabilities of the number of 
a, b, c-type edges, respectively. 

Here —a is the growth rate; a is called the surface tension, see below. 

What we've done above is a standard operation, called Legendre du- 
ality. Set a = 1 and set b = e x ,c = e Y . Then our expression for the 
normalized logarithm of the partition function is 

\ogZ = / log(l + e x z + e Y w) 



jui = i w i =1 2-Kiz2mw 
We can rewrite this integral as 

dz dw 



\ogZ = R(X,Y) = f j log(l + z + 

J \z\=e x J lui|=e y 



W) 



'\ z \=ex j\ w \ =e Y 27viz2mw 

Here R(X, Y) is called the Ronkin function of the polynomial P(z, w) = 
1 + z + w. 

We have shown in (jSJ) above that the surface tension is the Legendre 
dual of R: 

- a(s,t) = R(X,Y) - sX -tY, (6) 
where s = dR ^' Y "> and t = J^S ■ Both R and a are convex functions; R is 

dX dY 

defined for all (X, Y) G M 2 and o is only defined for (s, t) in the triangle of 
allowed slopes {(s, t) | < s,t,l — s — t < 1}. 
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Recall that 1 — s — t, s, t are proportional to the angles of the triangle 
with sides a,b,c. We have 

^— = X = \ogb/a 

as 

^ = Y = log c/a. (7) 

From this we can easily find 
Theorem 8. 

a(s, t) = — (L(tts) + L(wt) + L(7r(l - s - t))), 

7T 

where 

L(6) = - log 2 sin t eft 
Jo 

is the Lobachevsky function JW^ . 

Combined with ([6]) this gives an expression for R in terms of L as well. 
See Figures [[3] and [TH for plots of a and R. 



8 Boundary conditions 

As we can see in Figures H] and [SI boundary conditions have a large influence 
on the shape of a random dimer configuration. 

This influence is summed up in the following theorem. 

Theorem 9 (Limit shape theorem |5J). For each e > let 7 e be a closed 
curve in el? and which can be spanned by a stepped surface in el? . Suppose 
the 7 e converge as e — » to a closed curve 7. Then there is a surface S 
spanning 7 with the following property. For any 5 > 0, with probability 
tending to 1 as e — > a uniform random stepped surface spanning r y e lies 
within 5 of E . The surface S is the graph of the unique function h : 
Pm — > K. which minimizes 



min / o~(Vh)dx dy, 
h Ju 



w 
where 

cr(s, t) = (L(tts) + L(wt) + L(tt(1 - s - t))) 

7T 

and U is the region enclosed by the projection 0/7. 
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Figure 13: 
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Proof. Here is a sketch of the proof of this theorem. The space L(pf) of Lip- 
schitz functions with sloped in {(s,t) | < s, t, 1 — s — t < 1} having fixed 
boundary values (i.e. so that their graph spans 7) is compact in any natural 
topology, for example the L°° metric. The integral of a(Vh) is a lower semi- 
continuous functional on this space £(7), that is, given a convergent sequence 
of functions, the surface tension integral of the limit is less than or equal to 
the limit of the surface tension integrals. This follows from approximability 
of any Lipschitz function by piecewise linear Lipschitz functions on a fine 
triangular grid (see the next paragraph). The unicity statement follows from 
the convexity of o. 

Let L e (j) be the set of stepped surfaces spanning j e . Given a small 5 > 0, 
take a finite cover of £(7) by balls of radius 5. One can estimate the number 
of elements of L t {^f) contained in each ball as follows. This is essentially a 
large-deviation estimate. Given a function / G £(7), to estimate the number 
of stepped surfaces of L e (j) within 8 of /, we triangulate fl into equilateral 
triangles of side y/e. Because / is Lipschitz, Rademacher's theorem says that 
/ is differentiable almost everywhere. In particular on all but a small fraction 
of these mesoscopic triangles, / is nearly linear. The number of surfaces lying 
near / can then be broken up into the number of surfaces lying near / over 
each triangle, plus some errors since these surfaces must glue together along 
the boundaries of the triangles. It remains then to estimate the number of 
stepped surfaces lying close to a linear function on a triangle. This number 
can be shown, through a standard subadditivity argument, to depend only 
on the slope of the triangle and its area, in the sense that the number of 
stepped surfaces lying close to a linear function of slope (s, t) on a triangle 
of area A is exp(— Ao~(s, t)(l + o(l))) for some function o. 

It remains then to compute a(s, t), which is minus the growth rate of the 
stepped surfaces of slope (s,t). This was accomplished for the torus above, 
and again a standard argument shows that the value for the torus is the same 
as the value for the triangle or any other shape. □ 

Exercise 17 (up-right lattice paths). Consider the set of all up-right lattice 
paths in Z 2 starting at (0,0), that is, paths in which each step is (1,0) or 
(0, 1). Given a, b > consider a measure on up-right paths of length n which 
gives a path with h horizontal steps and v vertical steps a weight proportional 
to a h b v . What is the partition function for paths of total length n ? What is the 
typical slope of a path of length n for this measure? What is the exponential 

6 that is, the slope at almost every point 
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growth rate of unweighted paths with average slope s ? Describe the Legendre 
duality relation here. 

Exercise 18. On 7? , edges on every other vertical column have weight 2 
(other edges have weight 1). Redo the previous exercise if the paths are 
weighted according to the product of their edge weights times the a v b h fac- 
tor. 



9 Burgers equation 

The surface tension minimization problem of Theorem [9] above can be solved 
as follows. The Euler-Lagrange equation is 

div(Va(V/i)) = 0. (8) 

That is, any surface tension minimizer will satisfy this equation locally, at 
least where it is smooth. Here we should interpret this equation as follows. 
First, V/i = (s, t) is the slope, which is a function of x, y. Then a(Vh) defines 
the local surface tension as a function of x, y. Now Vcr is the gradient of o as 
a function of s,t. By Legendre duality, see ([7]), we have Vcr(s, t) = (X, Y). 
Finally the equation is that the divergence of this is zero, that is ^ + = 0. 
Substituting X = \ogb/a,Y = log c/a, see (|7j), this is 

dlog(6/q) + rflog(c/a) = Q , g . 
dx dy 

Consider the triangle with sides l,b/a,c/a, placed in the complex plane 
so that the edge of length 1 goes from to 1 as shown in figure [15j Define 




Figure 15: 

z, w G C to be the other two edges of the triangle so that 1 + z + w = and 
z = {b/a)e i ^- 9c) and w = (c/a)e l{7T+8b) . 
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Then 



V/i = (s, t) = (— , — ) = (- axg(-tu), - arg(--)). 

7T 7T 7T 7T Z 



In particular we have the consistency relation 
which gives 



Sy foxy foyx ^x-j 



I m f^ + ^ =0 . (10) 
V z w J 

Combining (jHJ) and ( fTOl) gives 

Theorem 10. For 2;, w as defined above, a general solution to the Euler- 
Lagrange equation is given by 

^ + ^ = 0. (11) 

z w 

This equation can be solved using the method of "complex characteris- 
tics". The solutions can be parametrized by analytic functions in two vari- 
ables. 

Corollary 11. For any solution to ffl\) there is a bivariate analytic function 
Qo for which 

Q (z,xz + yw) = 0. (12) 

Conversely, any equation of this type (with Qo analytic in both variables, 
which defines z, w implicitly as functions of x,y) gives a solution to 07]) . 

Proof. The existence of an analytic dependence between z and xz + yw is 
equivalent to the equation 

z x {xz + yw) x 



z y (xz + yw) y 

or 

z x (xz y + w + yuiy) = z y (xz x + z + yw x ). 
However since z and w are analytically related, z x w y = z y w x , leaving 

Finally, since z + w + 1 — 0, z y — —w y and so this last is equivalent to 
(HO. " □ 
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9.1 Volume constraint 

If we impose a volume constraint, that is, are interested in stepped surfaces 
with a fixed volume on one side, we can put a Lagrange multiplier in the 
minimization problem, choosing to minimize instead J J a + A J J h. This 
will have the effect of changing the Euler-Lagrange equation to 

div(V(r(V/i)) = c 
for a constant c. Equation ffTTj) then becomes 

^ + ^ = c (13) 

z w 



for the same constant c and equation (1121) becomes 

Q{e- cx z, e~ cy w) = 0. (14) 

For some reason the c ^ case has a more symmetric equation than the 
c = case. 



9.2 Frozen boundary 

As we have discussed, and seen in the simulations, the minimizers that we 
are looking for are not always analytic, in fact not even in general smooth. 
However Corollary [11] seems to give z, w analytically as functions of x, y. 
What is going on is that the equation (fTTj) is only valid when our triangle is 
well defined. When the triangle flattens out, that is, when z,w become real, 
typically one of the angles tends to it and the other two to zero. This implies 
that the slope (s, t) is tending to one of the corners of the triangle of allowed 
slopes, and we are entering a frozen phase. The probabilities of one of the 
three edge types is tending to one, and therefore we are on a facet. 

The boundary between the analytic part of the limit surface and the facet 
is the place where z, w become real. This is called the frozen boundary, 
and is described by the real locus of Q. 

Note that when z, w become real, the triangle can degenerate in one of 
three ways: the apex — w can fall to the left of 0, between and 1, or to 
the right of 1. These three possibilities correspond to the three different 
orientations of facets. 
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9.3 General solution 

For general boundary conditions finding the analytic function Q in ffTTl) which 
describes the limit shape is difficult. 

For boundary conditions resembling those in Figure HI however, one can 
give an explicit answer. Let Q be a polygon with 3n edges in the directions 
of the cube roots of 1, in cyclic order 1, e 27 ™/ 3 , e 47 ™/ 3 , 1, . . . , as in the regular 



hexagon or Figure [Tafl . Suppose that Q can be spanned by (a limit of) 




Figure 16: In this case Q has degree 3 and the frozen boundary is a cardiod. 

stepped surfaces. Suppose that there is no "taut edge" in fl, that is, every 
edge in the underlying graph has probability lying strictly between and 1 
of occurring in a dimer cover (an example of a region with a taut edge is the 
union of two regular hexagons, joined along a side). Then the limit shape 

7 The cyclic-order condition can be relaxed by allowing some edges to have zero length 
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arises from a rational plane curve 



Q(z,w)= 



C;jZ l W- 



0<i,j,i+j<n 



of degree n (or < n if there are edges of zero length). It can be determined 
by the condition that its dual curve (the frozen boundary) is tangent to the 
3n edges of Q, in order. 

Note that these polygonal boundary conditions can be used to approxi- 
mate any boundary curve. 

Exercise 19. Consider stepped surfaces bounding the regular hexagon. Show 
that plugging in Q (u,v) = 1 + u + u 2 — v 2 into fflfy gives the solution to the 
limit shape for this boundary. 

Let Q(u, v) = 1 + u 2 + v 2 + r(u + v + uv), where r > 2 is a parameter. 
Show that this value of Q in (T$) gives a solution to the volume- constrained 
limit shape, with volume which is a function of r. 

10 Amoebas and Harnack curves 

As one might expect, everything we have done can be generalized to other 
planar, periodic, bipartite graphs. Representative examples are the square 
grid and the square-octagon grid (Figure [9]). 

For simplicity, we're going to deal mainly with weighted honeycomb 
dimers. Since it is possible, after simple modifications, to embed any other 
periodic bipartite planar graph in a honeycomb graph (possibly increasing 
the size of the period), we're actually not losing any generality. We'll also 
illustrate our calculations in an example in section 110.41 

So let's start with the honeycomb with a periodic weight function v on 
the edges, periodic with period £ in directions x and y. As in section H~2l we 
are led to an expression for the partition function for the nl x ni torus: 




where 



Z. 



II II P{z,w), 



(-l)n w n =(-i) T 2 
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and where P(z, w) is a polynomial, the characteristic polynomial, with co- 
efficients depending on the weight function v. The polynomial P(z, w) is the 
determinant of K(z, w), the Kasteleyn matrix for the £ x £ torus (with appro- 
priate extra weights z and w on edges crossing fundamental domains); as such 
it is just the signed sum of matchings on the £ x £ torus consisting of a single 
£ x £ fundamental domain (with the appropriate weight (— \} hxhy z hx w hy ). 

The algebraic curve P(z, w) = is called the spectral curve of the 
dimer model, since it describes the spectrum of the K operator on the whole 
weighted honeycomb graph. 

Many of the physical properties of the dimer model are encoded in the 
polynomial P. 

Theorem 12. The free energy is 

^ , 1 i / t t \ 1 f .dzdw 
-\ogZ:= - km —\ogZ{H ne ) = -y— ^ / logP (z,w) . 

The space of allowed slopes for invariant measures is the Newton polygon of 
P (the convex hull in M 2 of the set {(i,j) : z l w l is a monomial of P}). The 
surface tension a(s, t) is the Legendre dual of the Ronkin function 



R(X,Y) = 7^[ I log P(z,w) 

™y j\ z \= e x ju 



dz dw 

(2vri) 2 J {zl=e x J H=eY ~*~ v ~'^ y z w ' 



We'll see more below. 



10.1 The amoeba of P 

Tke amoeba of an algebraic curve P(z, w) = is tbe set 

A(P) = {(log |4 log H) G M 2 : P{z,w) = 0}. 

In other words, it is a projection to M 2 of the zero set of P in C 2 , sending 
(z, w) to (log \z\, log Note that for each point (X, Y) e M 2 , the amoeba 
contains (X,Y) if and only if the torus {(z,w) G C 2 : \z\ = e x , = e Y } 
intersects P(z,w) = 0. 

The amoeba has "tentacles" which are regions where z — > 0, oo, or 
w — > 0, oo. Each tentacle is asymptotic to a line a log \z\ + f3 log \w\ + 7 = 0. 
These tentacles divide the complement of the amoeba into a certain number 



45 



of unbounded complementary components. There may be bounded comple- 
mentary components as well. 

The following facts are standard; see [221 [19]. The Ronkin function of P 
is convex in IR 2 , and linear on each component of the complement of A(P)H. 
The Legendre duality therefore maps each component of the complement of 
A(P) to a single point of the Newton polygon J\f(P). This is a point with 
integer coordinates. Unbounded complementary components correspond to 
integer points on the boundary of A/"(P); bounded complementary compo- 
nents correspond to integer points in the interior of 

See Figure [18] for an example of an amoeba of a spectral curve. 



10.2 Phases of EGMs 

Sheffield's theorem, Theorem O says that to every point (s, t) in the Newton 
polygon of P there is a unique ergodic Gibbs measure The local statistics 
for a measure /i Sj t are determined by the inverse Kasteleyn matrix K^y, 
where (X, Y) is related to (s, t) via the Legendre duality, S/R(X, Y) = (s, t). 
As discussed in section ISTTj values of K~ x are (linear combinations of) Fourier 
coefficients of 1/ P(z,w). In particular, if P(z,w) has no zeroes on the unit 
torus {\z\ = \w\ = 1}, then 1/P is analytic and so its Fourier coefficients 
decay exponentially fast. On the other hand if P(z,w) has simple zeroes on 
the unit torus, its Fourier coefficients decay linearly. 

This is exactly the condition which separates the different phases of the 
dimer model. If a slope (s, t) is chosen so that (X, Y) is in (the closure of) 
an unbounded component of the complement of the amoeba, then certain 
Fourier coefficients of 1/P (those contained in the appropriate dual cone) 
will vanish. This is enough to ensure that /j s j is in a frozen phase (yes, this 
requires some argument which we are not going to give here). For slopes 
(s,t) for which (X, Y) is in (the closure of) a bounded component of the 
complement of the amoeba, the edge-edge correlations decay exponentially 
fast (in all directions) which is enough to show that the height fluctuations 
have bounded variance, and we are in a gaseous (but not frozen, since the 
correlations are nonzero) phase. 

In the remaining case, (X, Y) is in the interior of the amoeba, and P 
has zeroes on a torus. It is a beautiful and deep fact that the spectral 

8 This shows that the complementary components are convex. 

9 Not every integer point in Af(P) may correspond to a complementary component of 

A. 
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curves arising in the dimer model are special in that P has either two zeros, 
both simple, or a single noda^l over each point in the interior of A(P). As 
a consequence^! in this case the edge-edge correlations decay quadratically 
(quadratically in generic directions — there may be directions where the decay 
is faster). It is not hard to show that this implies that the height variance 
between distant points is unbounded, and we are in a liquid phase. 

10.3 Harnack curves 

Plane curves with the property described above, that they have at most two 
zeros (both simple) or a single node on each torus \z\ = constant, \w\ = 
constant are called Harnack curves, or simple Harnack curves. They 
were studied classically by Harnack and more recently by Passare, Rullgard, 
Mikhalkin, and others [221 H5] . 

The simplest definition is that a Harnack curve is a curve P(z, w) = 
with the property that the map from the zero set to the amoeba A(P) is at 
most 2 to 1 over A(P). It will be 2 to 1 with a finite number of possible 
exceptions (the integer points of Af(P)) on which the map may be 1 to 1. 

Theorem 13 ([T6J [15]). The spectral curve of a dimer model is a Harnack 
curve. Conversely, every Harnack curve arises as the spectral curve of some 
periodic bipartite weighted dimer model. 

In [15] it was also shown, using dimer techniques, that the areas of com- 
plementary components of A(P) and distances between tentacles are global 
coordinates for the space of Harnack curves with a given Newton polygon. 

10.4 Example 

Let's work out a detailed example illustrating the above theory. We'll take 
dimers on the square grid with 3x2 fundamental domain (invariant under 
the lattice generated by (0,2) and (3,1)). Take fundamental domain with 
vertices labelled as in Figure [T7] — we chose those weights to give us enough 
parameters (5) to describe all possible gauge equivalence classes of weights on 
the 3x2 fundamental domain. Letting z be the eigenvalue of translation in 

10 A node is point where P = looks locally like the product of two lines, e.g. P(x, y) = 
x 2 - y 2 + 0(x, y) 3 near (0, 0). 

n We showed what happens in the case of a simple pole already. The case of a node is 
fairly hard. 
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Figure 17: 



direction (3, 1) and w be the eigenvalue of translation by (0, 2), the Kasteleyn 
matrix (white white vertices corresponding to rows and black to columns) is 



K 



1 + i 

w 

C 



w 



d 

-b+ 1 



We have 
P(z,w) 



det K(z, w) 

l + b + ab + bc + d + e 



1 + a + ab + c + d + ae a , ce , z 
+ — - bw + — + d—. 

w w z z w 



This can of course be obtained by just counting dimer covers of Z 2 /{(0, 2), (3, 1)} 
with these weights, and an appropriate factor (— l)^^ z t w^ when there are 
edges going across fundamental domains. Let's specialize to b = 2 and all 
other edges of weight 1. Then 



P(z,w) = 9-2w + 



z 



7 1 

— + - + — 

w z w 



The amoeba is shown in Figure [18j There are two gaseous components, 
corresponding to EGMs with slopes (0, 0) and (0, —1). The four frozen EGMs 
correspond to slopes (1, —1), (0, 1), (0, —2) and (—1, 0). All other slopes are 
liquid phases. 

For generic positive a, b, c, d, e there are two gas components. If we take 
c = 2,6 = |(3 ± y3) and all other weights 1 then there remains only one 
gaseous phase; the other gas "bubble" in the amoeba shrinks to a point and 
becomes a node in P = 0. There is a codimension-1 subvariety of weights 
for which P = has a node. Similarly there is a codimension-2 subvariety in 
which both gas bubbles disappear. For example if we take all edge weights 
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Figure 18: 



I then both gaseous phases disappear; we just have the uniform measure on 
dominos again. Dimer models with no gas bubbles, that is, in which P has 
genus zero, have a number of other special properties, see [TBI [To] . 

Exercise 20. How would one model the uniform square- grid dimer model 
inside the weighted honeycomb dimer model? 

Exercise 21. On the square grid with uniform weights, take a fundamental 
domain consisting of vertices (0, 0) and (1, 0). Show that (with an appropriate 
choice of coordinates) P(z,w) = l + z+w — zw. Sketch the amoeba of P. How 
many complementary components are there ? What frozen configurations do 
they correspond to? 

Exercise 22. Draw the amoebas for P(z,w) = a + z + l/z + w + l/w for 
various values of a > 4 (when \a\ > 4 this is a Harnack curve). Show that 
when a = 4, P has a node at z = w = — 1. 

II Fluctuations 

The study of the fluctuations of stepped surfaces away from their mean value 
is not completed at present. Here we'll discuss the case of the whole plane, 
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and stick to the case of uniform honeycomb dimers with the maximal mea- 
sure. Similar results were obtained |9] for uniform dimers on the square 
lattice, and (harder) uniform weights on these lattices for other EGMs [T2] . 
This last paper also computes the fluctuations for more general boundary 
conditions, in particular when there is a nontrivial limit shape. 



11.1 The Gaussian free field 

The result is that the fluctuations are described by a Gaussian free field. 
Here we discuss this somewhat mysterious object. 

The Gaussian free field in two dimensions is a natural generalization of 
one-dimensional Brownian motion. Recall that a Brownian bridge B(t) is a 
Brownian motion on R started at the origin and conditioned to come back 
to its starting point after time 1. It is a random continuous function on [0, 1] 
which can be defined by its Fourier series 

oo 

B(t) = V^sin(TrH), 
^— ' k 

k=l 

where the coefficients are i.i.d. standard normals. We can consider the 
Brownian bridge to be a Gaussian measure on the infinite-dimensional space 
L 2 ([0, 1]), and in the basis of Fourier series the coefficients are independent. 

The Gaussian free field on a rectangle has a very similar formula. We 
write 

oo 

GFF(x, y) = / . 2 J ' fc 2 sin(7rjx)sin(7r%), 

j,k=i V 3 + « 

where Cjk are i.i.d. standard normals. This defines the GFF as a Gaussian 
measure on a space of functions on [0, l] 2 . One problem with (or maybe 
we should say, feature of) this description is that the above series converges 
almost nowhere. And indeed, the GFF is not a random function but a random 
distribution. Given any smooth test function ip, we can define GFF(^) by 

oo „ 

GFF(^) = i o' k ; o / 4>(x,y)sm(njx)sm(Tiky)dxdy, 

~T, \1 + k J\o,i] 2 



and this sum converges almost surely as long as ip is smooth enough (in fact 
ip continuous suffices). Asking where the GFF lives, exactly, is a technical 
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question; for us it suffices to know that it is well defined when integrated 
against smooth test functions. 

The Gaussian free field on an arbitrary simply-connected bounded planar 
domain fl has a similar description: it is a Gaussian process on (distribu- 
tions on) Q with the property that, when expanded in the basis of orthonor- 
mal eigenfunctions of the Laplacian, has coefficients which are independent 
normals with mean zero and variance 1/|A|, where A is the corresponding 
eigenvalue. 

An alternative and maybe simpler description is that it is the Gaussian 
process with covariance kernel given by the Dirichlet Green's function g(x, y). 
That is, the GFF on Vt is the (unique) Gaussian measure which satisfies 

E(GFF(z!)GFF(z 2 )) =g(z 1 ,z 2 ). 

From this description we can see that the GFF is conformally invariant: given 
a conformal mapping 0: Q — > Q', the Green's function satisfies g(zi,z 2 ) = 
g((f>(zi), (f)(z 2 ) ). This is enough to show that 

/ GFF (z)^(z)\dz\ 2 = [ GFF(w)V>(0 -1 H)l<M 2 , 
with the equality holding in distribution. 



11.2 On the plane 



The GFF on the plane has a similar formulation, but it can only be integrated 
against L 1 functions of integral zero. We have 

E((GFF(*i) - GFF(z 2 ))(GFF(z 3 ) - GFF(z 4 ))) = 

= g(zi, Z3) - g(z u z 4 ) - g(z 2 , z 3 ) + g(z 2 , z 4 ) 



{zi ~ z 3 ){z 2 - Z A ) 



(zi - z 4 )(z 2 - z 3 ) 



where the Green's function g(z±, z 2 ) = log \ zi — z 2 \. 



11.3 Gaussians and moments 

Recall that for a mean-zero multidimensional (even infinite dimensional) 
Gaussian process X, if linear functions of X then E(xi • • -x n ) 
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is zero if n is odd and if n is even then 

E(xi ■ • -X n ) = ^E(x CT (i),X CT ( 2 )) • ..^{x a (n-i),X a ( n )), (15) 

where the sum is over all pairings of the indices. For example 

E(xi • • • x 4 ) — E(xi, x 2 )E(x 3 , x A ) + E(xi, a; 3 )E(a;2, x 4 ) + E(a?i, x 4 )E(x 2 , x 3 ). 

This shows that the moments of order two, ~E(xiXj), where Xi, Xj run over 
a basis for the vector space, determine a Gaussian process uniquely. Another 
fact we will use is that any probability measure whose moments converge to 
those of a Gaussian, converges itself to a Gaussian pQ. 

11.4 Height fluctuations on the plane 

We show here that the height fluctuations for the measure u = ui i on dimer 

3 ' 3 

covers of the honeycomb converge to the Gaussian free field. This is accom- 
plished by explicitly computing the moments E((h(zi) — hfo)) . . . (h(z n -i) — 
h(z n ))) and showing that they converge to the moments of the GFF. 

In fact we will only do the simplest case of the first nontrivial moment. 
The calculations for higher moments are similar but more bookkeeping work 
is needed. 

Let's fix four points z±, Z2, z%, Z4 G C and for each e > let v±, V2, V3, V4 be 
faces of eH, the honeycomb scaled by e, nearby. 

In this section for convenience we will use the "symmetric" height func- 
tion, where u is 1/3 on each edge. To compute h{v\) — h(v 2 ), we take a path 
in the dual graph from v\ to t>2 and count the number of dimers crossing it, 
with a sign depending on whether the dimer has white vertex on the left or 
right of the path. The height difference h(vi) — h(v2) is this signed number of 
dimers, minus the expected signed number of dimers. When e is small, v±,v 2 
are many lattice spacings apart and we can choose a path which is polygo- 
nal, with edges in the three lattice directions. By linearity of expectation it 
suffices to consider the case when both paths from v\ to v 2 and from v 3 to v 4 
are (disjoint) straight lines in lattice directions. So let us consider first the 
case when both lines are vertical. 

Let ai, . . . , a n be the edges crossing the first line (the line from v\ to u 2 ), 
and b\, . . . , b m be the edges crossing the second line (the line from v 3 to t> 4 ). 
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Then 



(Id ai -^(Id 6 ,-^) 



E[(h( Vl ) - h(v 2 ))(h(v 3 ) - h(v 4 ))\ = J2J2 E 

i=i j=i 

n m 

= EE E ( Id - Id O"q 



i=l j=l 

Here Id fli is the indicator function of the presence of edge aj, and E(Id a J = 
E(Id b:i ) = |. This moment is thus equal to 

n m 

i=i j=i 

where = w^bj and bj = w^b^. 

At this point we need to use our knowledge of i^ _1 (b,w) for points b, w 
far from each other. We have 

Lemma 14. 

2TTi(x-y)/3 



/ e 2ni(x-y)/3 \ 

K-\wo,o, Kv) = eRe [ ^-j^ - ^ j + 0{ 



\x\ + \y\ 



Here the e comes form the scaling of the lattice by e. 

This lemma is really the fundamental calculation in the whole theory, 
so it is worth understanding. First of all, recall that the values of K _1 (b,w) 
are just the Fourier coefficients for the function l/(z + w + 1). As mentioned 
earlier, if this function were smooth on the unit torus \z\ = \w\ = 1, its 
Fourier coefficients would decay rapidly. However 1 / (z+w+1) has two simple 
poles on the torus: at (z,w) = (e 2n ^ 3 , e 47 ™/ 3 ) and its complex conjugate 
(g47™/3^ e 27n/3)_ The Fourier coefficients still exist, since you can integrate a 
simple pole in two dimensions, but they decay only linearly. Moreover, for 
\x\ + \y\ large the (x, y)-Fourier coefficient only depends on the value of the 
function 1/P near its poles. Indeed the coefficient of the linearly-decaying 
term only depends on the first derivatives of P at its zeros. This implies 
that the large-scale behavior of K 1 — and hence the edge-correlations in the 
dimer model — only depend on these few parameters (the location of the zeros 
of P and its derivatives there). 
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With this lemma in hand we can compute (when the lattice is scaled by 

e) 

e 2 £ 2 e 2 e Ani(x-y)/3 £ 2 e -4ni(x-y)/3 

^ * ^ 4tT 2 (u 2 — Ml) 47T 2 (M 2 — Ml) 2 4:TT 2 \U2 — Ml| 2 47T 2 |w 2 — Mi | 2 ^ 

+0(e 3 /K- M2 | 3 ), 

where Mi is a point near a, and m 2 a point near fej. Summing over i,j, the 
terms with oscillating numerators are small, and this becomes 



21 J Z-i 



2Re / / -— ? -d Ul du 2 + 0(e) (16) 

47r z (Mi — m 2 ) 



1 teg ^"^"^ +0(6). 



27T 2 (^f2-^3)(^l-«4) 

Remarkably, we get the same integral ( fl~6j) when the paths are pointing in the 
other lattice directions, even when they are pointing in different directions. 

12 Open problems 

We have discussed many aspects of the dimer model. There are many more 
avenues of research possible. We list a few of our favorites here. 

1. Height mod k. What can be said about the random variable e iah where 
a is a constant and h is the height function? This is an analog of the 
spin-spin correlations in the Ising model and is a more delicate quantity 
to measure that the height function itself. Standard Toeplitz techniques 
allow one to evaluate it in lattice directions, and it is conjectured to 
be rotationally invariant (for the uniform square grid dimers, say). See 
[21] for some partial results. Can one describe the scaling limit of this 
field? 

2. All-order expansion. How accurately can one compute the partition 
function for dimers in a polygon, such as that in Figure 1161 .'' For a 
given polygon, the leading asymptotics (growth rate) as e — > is given 
by somewhat complicated integral (that we don't know how to evalu- 
ate explicitly, in fact). What about the asymptotic series in e of this 
partition function? For the random 3D Young diagram, this series is 
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important in string theory. (Note that the partition function for the 
uniform honeycomb dimer in a regular hexagon has an exact form, see 
Exercise |7]) . 

3. Bead model and Young tableaux. See [I]. For the a, b, c- weighted hon- 
eycomb dimer, consider the limit b = c = l,a — > 0. Under an appro- 
priate rescaling the limit is a continuous model, the bead model. The 
beads lie on parallel stands and between any two bead on one strand 
there is a bead on each of the adjacent strands. This model is closely 
related to Young tableaux. Can one carry the variational principle 
over to this setting, getting a limit shape theorem (and fluctuations) 
for random Young tableaux? 

4. Double-dimer model. Take two independent dimer covers of the grid, 
and superpose them. Configurations consist of loops and doubled edges. 
Conjecturally, in the scaling limit these loops are fractal and described 
by an SLE4 process. In particular their Hausdorff dimension is conjec- 
tured to be 3/2. 
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